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ABSTRACT 


This  report  contains  methods  for  constructing  iso-damage  curves  for 
structures.  The  results  for  the  cylindrical  shell  are  given  in  detail. 
It  is  shown  that  by  starting  with  the  theory  presented  here,  we  can  de¬ 
rive  the  empirical  relation  developed  by  Johnson  several  years  ago.  The 
theory  of  damage  due  to  short  duration  contact  explosion  is  presented 
and  the  results  are  compared  with  experiment.  A  series  of  curves  are 
presented  which  give  the  damage  sensitivity  of  cylinders  as  a  function 
of  the  ratios  of  diameter  to  thickness  and  length  to  diameter. 
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I.  Introduction  and  background 

The  iso-damage  concept  was  originated  at  Ballistic  Research  Laboratories 
in  the  early  1950's.  Just  as  isothermal  means  constant  temperature,  so 
iso-damage  means  constant  damage.  The  iso-damage  curve  consists  of  a 
plot  of  incident  pressure  as  ordinate  with  incident  impulse  as  abscissa 
such  as  shown  in  Fig.  1. 


i 


Fig.  1  Typical  Iso-Damage  Curve 


The  experimental  iso-damage  curve  is  formed  by  testing  a  structure  un¬ 
der  a  series  of  explosive  weights  and  distances  of  explosion  to  the  tar¬ 
get.  For  each  explosive  weight  and  distance  a  certain  damage  level  will 
occur.  A  curve  which  is  faired  through  the  same  damage  level  for  var¬ 
ious  impulse  and  pressure  is  called  an  iso-damage  curve.  Each  struc¬ 
ture  will  have  a  series  of  these  curves,  one  for  each  damage  level.  The 
asymptotes  shown  by  the  dotted  lines  are  minimum  values  of  pressure  and 
impulse  for  which  damage  level  A  will  occur.  Thus  for  -Z  <  _ro  ,  no 
matter  what  the  value  of  P  there  will  be  no  damage  at  "Damage  Level  A. " 
Likewise  for  P<  no  matter  what  the  level  of  I  there  will  be  no 
damage  at  "Damage  Level  A."  There  will  be  a  series  of  asymptotes  —  a 
pair  for  each  damage  level. 

II.  Iso-damage  theory 
A.  General  concepts 

In  a  previous  report‘d*  iso-damage  theory  was  illustrated  by  using  the 
exponential  decay  curve  as  a  first  approximation.  For  example,  assume 
that” the 'pressure  time  relation  is  as  follows: 

* Superscripts  refer  to  references  listed  at  the  end  of  the  report. 
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=  /'  €  '  "  [1] 

where  P  denotes  the  maximum  pressure  and  is  the  time  constant  of  the 
decay.  The  energy  flux  in  the  explosion  at  any  point  (i.e.  the  energy 
per  unit  area)  is  given  by  the  expression 

B,  =  -L  f dt  =  [2) 

where  fa  is  the  density  of  the  medium  and  C0  is  the  sound  velocity  in 
the  medium.  The  impulse  per  unit  area  is  given  by  2 

X  -  f  dt  =  P  ci 


Thus 


Ef  =  W 


Neglecting  any  dissipation  effects,  the  conservation  of  energy  demands 
that  the  energy  absorbed  by  the  structure  in  deforming  be  equal  to  the 
energy  directed  to  the  structure  from  the  explosion.  The  damage  that 
occurs  in  the  structure  can  be  measured  by  the  energy  absorbed  by  the 
structure  in  deforming.  Thus 

V  =  5  [5] 

where  V  is  the  total  energy  absorbed  by  the  structure  and  E  is  the 
energy  directed  to  the  structure  from  the  explosion.  The  total  energy 
£  can  be  written  in  terms  of  the  energy  flux,  Ef  by  the  expression 

E  -  Ef  A  [6] 

where  A  is  the  projected  area  of  the  structure  directed  toward  the  ex¬ 
plosion. 

The  energy  absorbed  by  the  structure  is  equal  to  the  internal  work 
done  by  the  structure  in  deforming.  For  a  given  level  of  damage  (i.e. 
a  given  deflection  distribution)  there  is  a  single  value  of  internal 
work  done  by  the  structure.  This  means  that  the  structure  does  a 
given  amount  of  internal  work  which  then  results  in  a  given  plastic  de¬ 
flection  distribution.  Thus  the  damage  in  the  structure  is  measured  by 
the  value  of  V  .  It  is  true  that  we  can  obtain  the  same  value  of 
for  different  deflection  distributions.  However  for  certain  types  of 
loads  and  structural  geometries  the  patterns  of  plastic  deflection  are 
fixed.  The  maximum  deflection  in  this  fixed  pattern  of  deformation  is 
a  measure  of  the  magnitude  of  V  and  therefore  of  damage  to  the  struc¬ 
ture.  Thus 

£>  =  V/A  [7] 

Substituting  into  [4] 

PI  =  2AcaV  I8! 

A 
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For  the  exponentially  decaying  pressure  the  iso-damage  curves  are  hyper¬ 
bolas  as  shown  in  Fig.  2: 


v 

A 


as  a  Parameter 

3 

B.  Conceptual  comparison  with  the  Johnson  Theory 

Recently  0.  T.  Johnson  of  BRL  offered  a  new  theory  of  blast  damage. 
Although  Johnson  derived  his  theory  from  empirical  considerations,  it 
can  be  shown  that  his  theory  is  fundamentally  well  grounded.  To  prove 
this  we  start  with  the  form  of  the  iso-damage  curve  derived  in  the  last 
section,  i.e.  _ 

PZ  =  A  Ef  [9] 

where  A  is  a  constant 

Consider  two  tests,  one  with  side  on  pressure  P,  ,  impulse  X,,  weight  W, 
and  distance  from  explosion  to  target  P,  ;  the  other  with  pressure  P. , 
impulse  Xz  ,  weight  and  distance  Pz  .  In  order  for  the  same  dam¬ 
age  to  occur  in  both  tests  the  same  energy  has  to  be  absorbed  in  the 
structure.  Thus,  for  the  same  damage 

P,  r,  =  P.  [io] 

4 


According  to  Cole  the  pressure  and  impulse  can  be  represented  as  a 
function  of  the  explosive  weight  and  distance  by  the  general  relations 

p*i(7?jy  _  ini 

x  =  2  l'y3)/S  [i2] 

where  ^  X  o7 }  ft  are  empirical  constants. 

Substituting  [11]  and  [12]  into  [10]  we  obtain  u 

—  /  w''1  /  t \  ft  —  s  '3  \  —  A  /  Wj.  3 ) 

■i(%  )j  KYx  )  =  -A.  ( V.  )  *  K  (  -Hi) 


S  -  /  — '  )  T?^r> 

'  (  ^4/ 


-  j- 


Thus 


[13] 

[14] 


From  Goodman's  curves  we  obtain  approximate  values  for  °<  and  i.e. 

£  &  /  }  /2  -x  / 

ThUS  \  ~  \'r 

(  (  'v'.-.)  [15] 

If  we  let  W  -  too  be  the  reference  weight  and  let  ^ -iv;  /<?,  /£ 

be  the  charge  weight  and  distance,  then 

(C‘~°//2.  )  ‘  ^  [16] 
or  (e'°°/sL  )  «  /OtA/~'S'  [17] 

3 

This  relation  is  very  close  to  the  one  given  by  Johnson.  It  just  var¬ 
ies  in  the  constant  and  the  power  of  .  If  we  had  used  more  accur¬ 
ate  values  for  ^  and  /£  we  would  have  obtained  values  even  closer  to 
Johnson's.  In  principle  it  is  seen  that  by  starting  with  the  energy 

concept  of  iso-damage  we  end  up  with  the  Johnson  Blast  Relationship.3 
* 

C.  More  accurate  formulation  of  iso-damage  curves  for  blasts  at  a  dis¬ 
tance  from  the  target 
6—8 

Brode  obtained  the  pressure  time  relation  for  spherical  blast  waves 
as  well  as  impulse  values  for  the  positive  phase  of  the  explosion.  As 
pointed  out  in  an  earlier  reference  the  maximum  underpressure  in  the 
negative  phase  of  the  pressure  is  generally  much  smaller  than  the  peak 
overpressure  at  the  shock  front  so  that  less  plastic  damage  will  occur 
in  the  negative  phase  than  in  the  positive  phase.  We  will  therefore 
use  the  positive  phase  impulse  and  energy  values  to  determine  the  dam¬ 
age.  Further  study  of  the  effect  of  the  negative  phase  on  the  plastic 
deformation  of  nonlinear  structures  will  certainly  be  warranted  at  a 
future  date. 

Although  Brode ' s  values  of  pressure  and  positive  duration  are  slightly 
in  error  (as  compared  with  experiment10)  for  spherical  pentolite  with  , 
values  of  greater  than  20  (/?  being  the  distance  of  target 

to  explosion  and  W  being  the  weight  of  explosive)  it  is  still  felt  that 
the  form  of  the  pressure  time  history  in  the  positive  phase  is  given 
accurately  by  Brode  and  therefore  his  impulse  values  will  be  satis¬ 
factory.  The  impulse  and  energy  that  will  be  used  to  describe  the  iso¬ 
damage  curves  are  the  positive  duration  values.  Thus 

T-X  fiRi 


=  jtc.  X 

where  T  is  the  positive  duration  of  the  overpressure, 


*  i.e.  More  accurate  than  the  theory:  -  Ref.  1 


Let 

?  =  i/T 

[20] 

then 

x  -  7  fa 

[21] 

and 

[22] 

The  pressure-time  relation  during  the  positive  phase  can  be  written  in  the 
form  6,  11  -<v  2 

PCs.  )  -  Ps  O  -  2:  )  <Z  [  23] 


where  is  a  nondimensional  constant  which  is  dependent  upon  /§  >  the 
peak  value  of  overpressure.  The  impulse  and  energy  can  then  be  written 


where 

tions. 


-r=  [24] 
Ef  =  ']£  L  ff'-Z)  [25] 

SO  X/r.T  -  A  C%)  =  fl  'r,-  2.)  [26] 

Eif°Co/r-p^  -  =  j^'rrz-z) e~* *"0^ *  i 27] 


•f,  ( Ps)  and  )  are  nondimensional  impulse  and  energy  func- 

Integrating,  we  obtain 

C Ps  )  -  ~  zfl  €  )  [28] 


X 

vi 


CPs)  =&[*-&  O  -  e-~)J 


[29] 


where  the  time  constant,  c*  ,  is  a  function  of  ‘s  •  Using  the  values 
of  impulse,  overpressure  and  positive  duration  given  by  Brode  '  ,  a  curve 
of  t,  C Ps)  as  a  function  of  Ps  has  been  obtained  and  is  given  in  Fig.  3. 

‘A  few  experimental  values  as  given  by  Goodman5  are  plotted  on  the  curve. 

For  low  overpressures  the  agreement  between  theory  and  experiment  is 
satisfactory.  However  for  the  higher  overpressures  (greater  than  10  at¬ 
mospheres)  there  are  too  few  experimental  points  and  too  much  scatter  to 
draw  any  definite  conclusions. 


The  value  of  ©C  was  determined  by  comparing  the  values  of  Ps  T  in 

Fig.  3  with  the  value  of  given  by  [28]  which  is  plotted  in  Fig.  4. 

The  nondimensional  energy  value,  C Ps  )  ,  was  then  obtained  from  this 

value  of  by  using  eq.  [29].  This  nondimensional  energy ' value  was  then 
plotted  in  Fig.  3. 


The  complete  iso-damage  curve  for  any  case  can  then  be  formulated  by  elim¬ 
inating  T  from  equations  [26]  and  [27],  i.e. 


Z  Pc£± 

Psfm  A(Ps)Ps*- 


[30] 
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CrPs)A^)4.rfi)  =  EfS.C' 


For  large  values  of  fj 


So,  for  large  Ps 


p  fps)  ^  £  r>- 
*  £  [>-*3 

*rPs}/fifPs) 


£PS 


2  fo  Co 


Equation  [34]  is  exactly  the  equatioi  obtained  for  the  simple  exponen¬ 
tial  in  Reference  1  and  in  Section  IIA  of  this  report. 

D.  Short  time  contact  explosions 

We  understand  blast  loading  as  that  loading  due  to  a  standoff  explosion 
which  produces  a  propagating  shock  wave  in  the  air.  There  are  other 
cases  of  impulse  response  under  explosive  loading  which  do  not  fall  into 
the  category  of  blast.  One  example  is  the  case  of  a  contact  explosion 
produced  by  sprayed  explosive.-*-2  The  sprayed  explosive  models  are  ex¬ 
posed  to  very  short  time  loading  of  various  distributions.  In  order  to 
derive  an  expression  for  the  response  in  this  case,  we  start  with  Ham¬ 
ilton's  Principle13 

^  ft,  ^  ~P  ~  O  [35] 

where  j  _  ^J^yO\(Cix+'ir*‘+s^r*‘)0lA,  the  kinetic  energy 

-  mass  per  unit  area  of  structure 

<AA  =  element  of  surface  area 

u,  a r si*s  =  velocities  in  the  three  coordinate  directions 

Assume  that  we  know  the  distribution  of  the  displacements  from  some  ex¬ 
perimental  work,  then  tA  ^UT  can  be  written 

U  -  (Aa  (t)Jz(A)  [36] 

AT  =  \T0  a  )  c A  ) 

-w  ~ **rc(-b)  f^r(A) 

where  f^r  are  the  distributions  of  6^  AaT  over 

the  surface  of  the  structure 

r  13 

The  variation  6  is  taken  exactly  as  in  the  elastic  problem.  This 

problem  is  equivalent  to  the  following  problem  in  the  Calculus  of  Var¬ 
iations: 

Find  the  functions  (x.)  j  ....  which 

take  on  given  values  for  sc  =  <x  and  sc  =  &  and  which  mini¬ 
mize  the  definite  integral 
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T~ 


yf_  F'  ■■■  ....  [37] 


The  result  is  that  F  must  satisfy  the  set  of  Euler  Equations 
In  our  case  F  -  7  ~  U 

?■  r  ^  Z'S'oj  #,'=  ^ 

Thus 

-  4  f s<*  ^ +  ~)  cd  A  —  V  '  +  UY  ' 

z  JA 


[38] 

[39] 

[40] 


where 


r  s 

V  =  J  (  J  o~.<dc:  )  d\/  (work  done  by  internal  forces)  [41] 
</•  ° 

[42] 


lA/  '  -  f  C X  U  f  )V  +  2.  yi^r  )  d  A 

-VI 

Thus  the  governing  equations  for  the  unknowns  ue  yj- 

””  j  C  ^  C 


are 


J^m  F  Va)  c/a  +  j^y'  -  £  xaaj  n  Ppa)  c/a 
'h  fyu  fy/CA)  dA  4  £j/'  -  f  YfAj/r)  f^rCA) 

A  JstT-  Ja 


d  A 


[43] 


f/*  ’U'CAie/A  +  £X.'  =  f  ZrA,-/)  jC  (A  )d A 

A  °-w~o  JA 

In  the  large  deflection  region  V^is  a  function  of  powers  of  \T0  *A~0 

so  that  the  functions  QZ'  <kV_'  <$_v' 

*  ow. 


)  0  iT. 


nonlinear  functions  of  oia  vc  /t^re,  •  The  equations  [43 j  are 

therefore  ordinary  nonlinear  differential  equations  for 

These  equations  can  be  written  in  the  standard  single  degree  of  freedom 


form 


14 


M  *•  +  &  X-  -  P('-k)  [44] 

where  M  is  the  generalized  mass,  Z.  is  the  generalized  resistance  and  P 
is  the  generalized  force.  Assume  that  (a  *J~  are  small  compared  to  -w 


For  the  lateral  deflection  «y*j- 


M  =  f^A)clA 

O  ur0 

?/-t)  ■=  /  Z/'A'i)fL-(A  )dA 


[45] 
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For  very  short  time  loading  eq.  [44]  can  be  simplified  considerably. 

If  a  dynamic  loading  P(-f)  is  applied  to  a  dynamic  system  with  one 
degree  of  freedom  (i.e.  a  system  satisfying  eq.  [ 44 [  the  external  work 
done  up  to  any  time  t  is  given  by 


-9- 


-  r&  y 

lA/'f-t  )  -  J  P  ft)  =  J*  pfi) 

The  velocity  x.  r-t)  is  determined  by  integration  of  eq.  [44] 

*(*)  -  pT  f  f  Fr+)  -  pro] 


Thus 


w'n)  -  £  PC-t)f\ fa[ p{±)-&r*S]J*J 


If  the  time  it  takes  to  reach  a  maximum  deflection  is  greater  than  the 
duration  of  the  load,  T,  then  eq.  [48]  need  only  be  integrated  up  to  T. 
Under  these  circumstances  Pr*.)  is  small  during  the  application 
of  the  load  and  can  therefore  be  neglected.  The  work  W  therefore 
becomes: 

m/'  -  f  pa)  £  -£■  f  prt)  Jt]  t49^ 


'  -  H2' 


iV  =  JZ.  1 

■2  A7 

_  r_ 

where  H  -  f  PftliJt 

H  is  the  total  impulse  of  the  external  load.  This  work,  w'done  by  the 
external  load  is  equal  to  the  work  done  by  the  internal  forces,  v/  in 
deforming  the  structure 

i*e*  v  -  R2 

2  F)  [5: 


V  - 


where 


V  =  f  ( 

v‘  ° 

H  ~  f  f  Z<'AJ£)£r('A)]ic/A 

A 


n  fAyo ,  f^yrA)]A 

Calculation  of  structural  energy  absorption 

The  loading  has  been  characterized  in  the  previous  sections  by  the  peak 
overpressure,  Ps  t  the  positive  duration  of  th^  overpressure,  Tt  and  the 
impulse  I.  For  either  blast  loading^-®'  4^  or  impulse  loading40,4® 
these  parameters  are  measurable.  The  main  question  is  concerned  with 
what  effect  the  magnitudes  of  overpressure,  duration  and  impulse  have 
on  the  damage  induced  in  the  structure.  The  characteristics  of  the 
structure  were  given  in  the  previous  sections  simply  by  V,  the  energy 
absorbed  in  deforming  the  structure. 

Suppose  there  are  direct  stresses  0"^  (p  G”i  and  shear  stresses 

acting  on  an  elemental  volume  of  the  structure  as  shown  in 

Fig.  5: 
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Pig.  Free  Body  Diagram  of  Elemental  Volume 

The  direct  and  shear  strains  are^  6  }  Y^  ,  YXi^  *  T^e 

stress  strain  curve  for  one  dimensional  strain  in  the  x.  direction  will 
be  as  shown  in  Fig.  6: 


Fig.  6  One  Dimensional  Stress-Strain  Curve 
The  elemental  work  for  straining  in  the  £  direction  will  then  be 

Jw  -  OZde^  t53l 

The  total  work  in  direct  stress  and  shear  will  then  be 


(d  W  -  (T^  d£x  +  (7^  dZy  +  ed  6  i 


+■  ’Z’x  +  tdYxi: 


[54] 


The  work  per  unit  volume  can  be  divided  into  the  work  done  in  changing 
the  shape  of  the  body  and  the  work  done  in  changing  the  volume  of  the 
body,  x.e. 


* '  -  3  J  ov<r  M7  r,  */£ 

f  z  ■> 


in  which 


CT*  =  ^  ^  +  <Ta- 

J 


=  «^<>c  jjLfg,  £*■ 


and  h 


are  the  octahedral  shear  stress  and  strain  given  by* 


^-o  -  j  \^('<rx.~  ~  )  t  (^i. ~(jZ)  +  6  r~  )  [  57  ] 

A  =  j  \j&e^*z?+&£3rclexf'H‘te*  -*teO' V 4  (eUr^+J^+efr^) 

If  the  material  is  incompressible  then  the  volume  does  not  change  and 
the  total  work  is 

|  [58] 

K 

For  an  incompressible  material  in  a  biaxial  stress  state  we  replace  the 
octahedral  shear  stress  and  strain  by  their  two  dimensional  counterparts 

91 ^ 

°I  and  e.i  .  The  work  done  per  unit  volume  in  distortion  is  then  '  ' 


"  =  /  07  J<c 


where 


<jl  -  v^^-cr;^  -t-<r7 *•+ 3  *• 


€■  =  2=. 


{j  V€**‘¥  ** 6 > y  ^  *  V 


The  total  work  done  through  the  entire  volume  is  then 

V  -  w  'dv'  [61] 

The  stress-strain  law  of  the  material  is  given  by  the  functional  rela¬ 
tionship  between  and  cfj  as  shown  in  Fig.  7 

I  ./<57=/Ve-) 


Fig.  7  Form  of  Stress  Strain  Law  of  Material 


* 


* 


The  stress  strain  law  can  have  a  variety  of  forms.  Most  practical  cases 
will  fall  under  the  following  categories: 

1.  Rigid-Linear  Hardening 

crL  =  <rs  +  k  ec 

K  —  s/apc 

X/  ATr  o  (pe-fce  +  1'/  p/ashc) 


Fig.  S  Rigid-Linear  Hardening  Law 

2.  Elastic-Linear  Hardening 


Ee-  c{<es 

07 -h  k  ce;  - es)  for  &■  >  <?s 


Fig.  9  Elastic-Linear  Hardening  Law 


3.  Elastic-Plastic  Power  Law 


Fig.  10  Elastic-Plastic  Power  Law 
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Using  each  of  these  laws,  substituting  into  [50]  and  integrating,  the  total 
energy  can  be  put  into  the  following  general  form 


where 

1.  Rigid-Linear  Hardening 

4'=  2 
£ / 
t'  --o 

2.  Elastic-Linear  Hardening 

£'=  / 

T'-O 


eL  *  +C3  c/')c/v  +  C+ 


c,  -  Vz 

4  ~  £e>r-  per  p  /e, 

4  =  ° 

4  =  ° 

c,  =  K/z. 
cx  =  (j^-/ces 
c3  =o 

4  =  *2#  aAzS 


[62] 


//e  frro 


where  A  =  surface  area  of  the  shell 
A.  =  shell  thickness 

3.  Elastic-Plastic  Power  Law 


4.  Bell  Law 


Q  -  p  +  / 

£'  -  O 

T '  -  p  +  i 

c  = 

•  />+  / 

Ci.  -  o 

<4=1  aA.cS 

a W 

/Z‘  -  o 

T'-O 

‘,-iA'O-t) 

Cz  -  o 

4 -o 

For  very  large  plastic  deformations  the  elastic  strains  are  very  small 
compared  to  the  plastic  strains  and  usually  4  can  be  placed  equal  to  0. 
The  energy  integral  [62]  is  given  in  terms  of  Cl  ,  which  is  a  function 
of  •  The  next  sections  of  the  report  will  be  devoted  to 

the  calculation  of  this  energy  for  particular  structures. 
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III. 


Energy  absorption  curves  for  various  failure  shapes  of  cylindrical 
shells 

A-  General  relations  for  energy  absorption  of  cylindrical  shells 


The  strains  £. 


of  the  biaxial  stress  field  can  be  written 


25 


3c  =■  2  tc, 


-  €u  -  Z  /rz 


Vxy.  -  y  -  2  Z  T 


[63] 

where  £  >  are  the  midsurface  strains,  2  is  the  radial  dis^ 

tance  from  the  midsurface  to  any  element  as  shown  in  Fig.  12  and 
K.  K,  '£■  are  the  curvatures  and  twist- 


are 
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£  -  du  +  L  Yrh. cfV 
'  “  ~JZ  Ok.  J 


£  -  J-  £^T_ 

1  G )cf 


■y~r 

(X. 


j.  J. 

+  Z  J 


[64] 


Y  -  +  JL  Old  +.  c >w  Jw 

U*.  <i.  Oef  Ok.  *Ocf 


The  curvatures  and  twist  are 
AT  - 
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0<t: 


4  = 


—  J- 


+  J~.^£ 

^  0c/> 


[65] 


Jc. 

For  very  large  deformations  under  intense  lateral  loading  the  mid¬ 
surface  strain  involving  (i.e.  the  nonlinear  terms)  should  pro¬ 
bably  be  greater  than  the  linear  terms  involving  and  .  Assum¬ 
ing  that  u  and  nT  and  their  derivatives  are  much  smaller  than  kkj~ 
and  its  derivatives,  we  have 


✓  _  J.  /  e)w  \  j-  <D5*r 
^x.  -  2  (  OkJ  ~  2  “ 


£=  r  1  )t±r_  z  Q^r 

<f  Ol.  tfi-  £  cp‘ 


-u/~ 
v 


[66] 


^  ..  2  2 

J"  c 0  X. 
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A  large  number  of  practical  cases  can  be  fitted  to  the  elastic  linear 
hardening  law  shown  in  Fig.  9.  Figure  13  shows  how  an  elastic-linear 
hardening  curve  can  be  fitted  to  an  actual  stress-strain  curve. 


Fig.  13  Fit  of  Elastic-Linear  Hardening  Law 

Temperature  can2glay  an  important  part  in  the  form  of  the  stress  strain 
law.  Figure  14  shows  the  effect  of  temperature  on  the  yield  stress 
and  hardening  characteristics  of  a  typical  elastic-linear  hardening  curve 


Fig.  14  Effect  of  Temperature  on  Stress-Strain  Law 

The  effect  of  increasing  temper' ture  is  to  decrease  the  yield  stress 
and  decrease  the  hardening. 


Since  the  elastic-linear  hardening  curve  can  be  used  to  describe  most  of 
the  critical  characteristics  of  the  metal,  we  will  limit  the  analysis  to 
this  type  of  stress-strain  curve.  Under  these  circumstances  the  energy 

a  "5  *\  mil  >  «  1  1  _ _  '  _  _ _ *  •  • _ “  y 


absorbed  in  a  pylindrical  shell  can  be  written 

J?  sltT 


S*  f  "a. 


■h  FA 


V7 


,e‘£3te> 


’  i  ^  '■ 


a  d  cf  a/ 


[6 


* 


where  A  -/-  K/e 

Substituting  the  expressions  for  the  strains  given  by  [66],  letting  x/jL 
and  ^ii-  =  ^w0  ■ffz.'jcp)  and  integrating  [67],  the  energy  absor¬ 


bed  can  be  written  in  the  following  convenient  form:^ 


The  integrals  are  dimensionless  quantities  which  are  functions  of  the 
dimensionless  ratios 

The  parameter  A  and  the  yield  stress  C7]T  are  outside  the  integrals. 
Therefore  the  value  of  the  integral  is  independent  of  both  the  hardening 
and  the  yield  stress. 
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There  are  three  main  types  of  failure  patterns  which  have  j^en  determin¬ 
ed  experimentally.  These  typical  deformation  patterns  '  are  shown 
in  Figures  15-19.  The  first  two  (Figs.  15,  16)  are  typical  patterns 
associated  with  blast  waves  originating  at  ci  distance  from  the  shell. 

The  other  figures  (i.e.  Figs.  17-19)  are  typical  patterns  from  contact 
explosions^ * -8  in  the  form  of  a  sheet  of  explosive  on  the  surface  of 
the  shell.  The  analytical  functions  describing  each  of  these  types  of 
patterns  and  the  corresponding  energy  absorbed  will  be  treated  in  the 
next  several  sections  of  the  report. 

B.  The  single  diamond  pattern  and  the  lobar  buckling  pattern 

1.  Criterion  for  determining  the  pattern 

31 

In  previous  work  a  criterion  was  determined  to  establish  whether 

the  hinged  single  diamond  pattern  (Fig.  15)  or  the  lobar  buckling 

pattern  (Fig.  16)  would  occur.  Over  the  past  several  years  it  has 

been  found  that  this  criterion  is  not  quite  accurate  enough.  Using 

the  same  ideas  as  in  the  previous  work,  i'c  will  be  assumed  that  the 

diamond  pattern  is  a  collapse  mode  and  that  the  shell  will  assume 

this  mode  of  failure  if  the  yield  condition  is  reached.  If  the 

load  which  produces  buckling  is  less  than  this  yield  load,  then 

buckling  will  occur.  In  order  to  calculate  the  yield  or  collapse 

load,  it  will  be  assumed  that  the  shell  is  thin  enough  to  take  the 

total  load  by  membrane  action  alone.  Under  these  circumstances  the 
,  ■  37 

stress  distribution  is  easily  found. 


pf'cf') 


Fig.  20  Statically  Loaded  Cylindrical  Membrane 

Let  the  length  of  the  shell  be  and  assume  that  it  is  supported 
by  a  diaphragm  at  each  end  (  21  -  iL  )  ,  if  the  origin  is  at 
the  center  of  the  shell  the  boundary  conditons  are 


~  O  /  X.  -  — 
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Fig.  15  Collapse  Pattern 
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PROFILE  OF  TEST  CYLINDER 


TEST  CYLINDER 

6061 -T6  ALUMINUM  ALLOY,  0.035"  WALL  THICKNESS 
3.00"  O.D.,  COSINE  SHEET  EXPLOSIVE  LOADING 


*1  Lonqitudinal  Distribution  Patterns 


[70] 


As/ /*/)  — 


A/xe/l  ~  -  'XL 

rrd»  -~  ^ 

^  ^  <?V 

At  the  same  time  of  contact  of  the  shell  with  the  shock  wave,  the  press¬ 
ure  will  be  a  maximum  on  the  side  of  the  shell  facing  the  explosion  and 
decrease  to  a  small  value  on  the  back  side.  The  load  distribution 
will  therefore  be  taken  in  the  following  form: 

Therefore  ‘  *  f 

a/**  -  ^ 

a/«.  = 

AJcf>  -  A  ft  g 

A  hinge  will  form  at  7c  zcfi-O  when  the  yield  condition  is  satisfied  at 
that  point,  i.e. 

i  A/^  1  -  Cf2'-/2'  [73] 

The  front  face  pressure  at  which  yield  (or  collapse)  will  start  is 
therefore: 

'p  -  _ L _ _  [74] 

/c  ^  - ~~~TT - — £ - - - 

v  C  Fa  *■  *  )  -t-  •*'  ^  1 

The  elastic  buckling  load  for  uniform  loading  of  the  cylinder  is  given 
by  Reynolds33 


where 


f6  ■=  Z*l£Jk- /£f  { 

3  Cp  —  0  L)  S  ^  y 


i/a^L 


-  A  2  3 


[ /#■£. 
A 


3  ~  2cf  ) 


4  X 


fee  -  $t  fe  '****£•  &<*  t  8 &  \ 

(  -  c*.  8  c  / 

^  '  ~  fe  /  <Vj_  -h  aA.  £«✓,_  \  „ 

■2-  l  t)e ~TZr7T77  J  >  -  /  . 


-  fe  /  fe/^\  I  J 

^  -  Ctn.  &c /l  )  '  fe  ~  J 4 

in  which  A_  is  the  cross  sectional  area  of  the  rings  which  supported  the 
cylinder  at  (see  Fig.  20),  b  is  the  width  of  the 


r 

t 


i 


* 


frame  in  contact  with  the  shell  (see  Fig.  21) 

j,  jL  *-•  h  H 


Fig.  21  Cross  Section  in  the  Vicinity  of  the  Ring  Support 
34 

Almroth  gives  correction  curves  for  buckling  under  nonuniform  loading. 
The  nonuniform  peak  buckling  pressure  can  be  written  as  '  where 

?e'*  *  fa  ,  *  >/  [17] 

2.  Comparison  of  the  collapse  and  buckling  relations  with  experiment 

A  large  nu^er  of  tests  were  carried  out  on  shells  of  various  sizes  by 
Schuman.  The  blasts  were  mostly  side  on  so  that  the  pressure  distri¬ 

bution  was  nonuniform  over  the  shell  circumference  but  uniform  over  the 
length.  Based  upon  previous  work35,36  an  assumption  of  ex’1-  /  for. 
calculation  of  the  collapse  load  seemed  reasonable.  The  buckling  load 
was  corrected  by  use  of  Almroth 1 s  curves.  Almroth  calculated  the  buck¬ 
ling  pressure  of  cylinders  with  the  following  pressure  distribution: 


7?  -  fo  +  <^-cf 


[78] 


Letting  c>  ‘  —  1°' 

f>o 

f,  - 

>  -f 

The  peak  load  is  at  c/>  -  o  and  has  a  value  of  { y  -p,  ) 

For  our  calculations  a  value  of  J>  ~  ,  7J~  was  used  to  estimate  the 
correction  due  to  nonuniform  loading. 

The  collapse  load  is  then  given  by 


/ 

y  / 

and  the  buckling  load  by 

fc'd. 

' 

K  £ 

3- 

[79] 


[80] 
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where  ^ /c 7~a  is  the  ratio  of  elastic  modulus  ^  yield  stress  and 

for  the  test  shells  is  given  in  a  previous  reference.  A;  is  the  Alm- 
roth  correction34  factor  for  nonuniform  loading.  Table  1  gives  the  re¬ 
sults  of  the  tneory  and  experiment.  The  last  two  columns  indicate  whether 
the  shell  collapsed  or  buckled  in  the  tests  and  what  was  predicted  from 


the  theory.  In  the  theory 

T-  HfS 

tra  &  <V  A. 

Collapse  will  occur 

T-f  :  >  7>a 

Buckling  will  occur 

In  Table  1  “ 

.  £>/t  -  *V2L 


Table  1  Collapse  and 

Buckling  Parameters 

B- Buckling 

ft  *- 

f&Jt 
or.  A- 

C-Collapse 

L'P 

Material  /<7a 

/•c 

<Te4- 

Theor . 

Exper . 

2 

158 

1040  Steel  1000 

1.25 

.38 

.44 

C 

C 

2.87 

158 

1.25 

.21 

.31 

C 

C 

3 

158 

1.25 

.20 

.29 

C 

C 

3.87 

158 

1.30 

.12 

.23 

C 

C 

4.87 

158 

1.35 

.08 

.19 

C 

c 

6 

158 

1.40 

.05 

.16 

C 

c 

8 

158 

1.50 

.03 

.13 

C 

c 

2.87 

86 

1.55 

.21 

.95 

C. 

c 

3 

86 

1.55 

.19 

.91 

C 

c 

6 

86 

1.60 

.05 

.47 

C 

c 

3 

316 

1.50 

.19 

.12 

B 

not  clear 

2.91 

172 

1.40 

.21 

.30 

C 

C 

3 

172 

1.40 

.20 

.29 

C 

c 

2.91 

79 

1.50 

.21 

1.03 

C 

c 

2.94 

158 

1.25 

.20 

.30 

C 

C 

2.94 

88 

1.40 

.20 

.81 

C 

C 

1.98 

176 

1.25 

.38 

.38 

C 

not  clear 

2.0 

1000 

Al  5052-48  333 

1.15 

.19 

.006 

B 

B 

3.0 

1000 

1.15 

.38 

.008 

B 

B 

5.0 

1000 

1.20 

.077 

.003 

B 

B 

3.0 

500 

1.25 

.20 

.017 

B 

B 

5 

500 

1.30 

.077 

.r  li 

B 

B 

7.67 

500 

1.35 

.033 

.007 

B 

B 

10 

500 

1.40 

.020 

.006 

B 

B 

3 

250 

1.25 

.197 

.005 

B 

B 

3 

125 

1.25 

.197 

.139 

B 

B 

-26- 


It  is  seen  that  with  very  few  exceptions  the  formulas  presented  in  the 
previous  section  accurately  predict  whether  the  shell  will  collapse  or 
buckle . 

C.  Energy  in  the  diamond  collapse  pattern 

The  plastic  energy  absorbed  in  the  collapse  pattern  during  large  defor¬ 
mation  was  obtained  in  an  earlier  report.3'1'  More  complete  curves  are  pre¬ 
sented  here  in  Figure  22  as  a  function  of  the  deflection  ratio 
which  is  a  more  convenient  parameter  than  the  one  used  for  collapse  in 
the  earlier  work.^  For  the  linear  hardening  curve  shown  in  Fig.  9  (see 
equation  67  for  X  ) ,  the  total  energy  can  be  written  as 

V  -  0~s  f  €LLL-±)  r,  +  x  v,  ~  A(/T rrcs 

~UT  i 

In  Fig.  22  the  nondimensional  energy  functions  x,  and  v,  are  plotted  for 
a  large  range  of  physical  parameters.  Note  that  these  functions  are  in¬ 
dependent  of  P/f  ana  for  L/p  >  g*  they  are  independent  of  u/o  .  Collap¬ 
se  involving  temperature  dependence  (see  Fig.  14)  and  other  hardening 
problems  can  be  completely  solved  by  this  set  of  curves  since  both  func¬ 
tions  X,  and  v,  are  independent  of  the  hardening  parameter  ^ 
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D.  Energy  in  the  buckled  pattern 

A  typical  buckling  pattern  that  could  occur  in  the  shell  due  to  side  on 
blart  is  shown  in  Figure  16.  The  theory  predicts  one  half  wave  length 
along  the  length  of  the  shell  and  waves  around  the  periphery,  where 

'vw  j.s  dependent  on  the  geometry  of  the  shell.  The  manner  of  buckling 
is  well  confirmed  by  experiment.  The  displacement  pattern  in  the  post 
failure  or  plastic  region  can  be  described  by  the  following  equation: 

-iZcfi 

slyj~  ~  'H  Cfi 

33 

The  number  of  peripheral  lobes,  /n<  is  given  by  Reynold's  as 

~  vAt hamt  ~  1841 

The  factor  of  1.23  is  not  exactly  correct.  Since  the  number  of  full 
waves,  'yv,  must  be  a  whole  number  this  factor  of  1.23  is  adjusted  so 
that  syo  is  the  whole  number  nearest  to  the  value  calculated  by  using 
the  factor  1.23.  The  parameter  -A.  varies  and  it  is  difficult  to  as¬ 
sign  a  reliable  number  to  this  parameter.  Nevertheless  a  range  for 
can  be  determined.  Based  upon  examination  of  experimental  result,  it 
seems  that  o<  A.<.  /  .  An  extensive  set  of  energy  absorption  curves 

for  A-  =•.  2-jr  is  shown  in  Figure  23-25.  These  curves  are  based  upon 
the  addition  of  membrane  and  bending  energy.  It  was  found  that  the 
membrane  was  much  greated  than  the  bending  for  buckling  in  this  large 
deflection  region.  Figures  26,  27  show  a  comparison  of  membrane  ab¬ 
sorption  energies  for  various  values  of  A.  between  0  and  1  for  several 
shell  geometries.  For  side  on  blasts  a  good  average  value  tor -A.  is  25. 

E.  Energy  in  the  short  duration  contact  explosion  pattern 


Some  typical  deformation  patterns  for  sprayed  explosive  loading  are 
shown  in  Figures  17-19.  The  deformation  patterns  vary  considerably, 
but  all  of  them  can  be  described  by  the  general  form 

^0[o e'Ccf O-** )  ^0] ;  «'=*/*  [85! 

The  constants  which  vary  with  each  pattern  are  G,  H,  C,  B,  K.  Figures 
28,  and  29  show  the  longitudinal  and  peripheral  displacement  distribu¬ 
tions  for  a  range  of  values  of  these  constant'..  Since  the  patterns  are 
always  symmetrical  about  &-o  and  7c  '  s ,  s'  onl^  one  half  of  the  pattern 
is  shown.  The  patterns  in  the  Stanford  Research  Institute  tests  vary 
greatly  from  those  in  the  Southwest  Research  Institute  tests  (see  Fig. 
19) .  It  is  therefore  difficult  to  conclude  what  a  typical  pattern  should 
be  for  a  given  value  of  Vjo  and  .  It  is  believed  that  these  short 

duration  explosive  tests  are  not  as  conclusive  as  the  collapse  and  buck¬ 
ling  phenomena  fi.om  stand  off  blast  as  presented  earlier  in  the  report. 
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Fig. 25  Energy  Absorption  in  the  Buckled  Pattern 
(k=. 25) 


It  was  found  that  the  energy  of  plastic  deformation  was  critically  de¬ 
pendent  upon  the  values  of  the  shape  parameters  G,H,C,B,K.  Therefore 
it  would  be  misleading  at  this  time  to  tie  down  energy  values  as  a 
function  of  l/q  D/±  .  However  in  the  next  section  of  the  report 

the  actual  impulse  values  obtained  from  tests  will  be  compared  with 
those  predicted  from  theory  by  using  representative  values  of  G,H,C,B,K. 

IV.  Impulse  response  of  cylindrical  shells  for  short  time  contact  explosions 

In  Section  IIB  of  this  report  there  is  a  derivation  of  the  impulse  as  a 
function  of  the  energy  absorbed.  For  the  cylindrical  shell  this  rela¬ 
tion  becomes  _ 

\  J  —  u  y  =  , 


where 


\/  -  IV  =  SL 
2./°? 

W  -  V  =■  work  done  by  internal  forces  (see  eq.  [68] 

-  - 

1  J.  L  *:'*,*» 


where  J CT  S*-,  &  )  is  the  impulse  distribution  and  X  is  the  peak  value 
of  impulse  per  unit  area.  Substituting  the  value  of  \/  from  [68] ,  H  and  M 
from  [87] ,  we  obtain  the  following  relation  for  a  rigid  —  perfectly 
plastic  material: 

( }  )  ~  f  >  yC> ;  )  [88] 

where  /  is  a  dimensionless  function  of  the  dimensionless  parameters 

is  the  maximum  deflection,  which  in 

all  cases  considered  here,  is  at  cf  =  o  3c.  -  .  Equation  [88] 

was  programmed  on  a  time  sharing  system  and  calculations  were  run  using 
a  cosine  distribution  of  impulse  over  180°  ^  using  representative  val¬ 
ues  of  the  constants  describing  the  deformation  pattern.  A  rigid  per¬ 
fectly  plastic  material  (see  Fig.  8)  was  assumed.  It  was  found  that  the 
following  values  of  the  constants  gave  a  deflection  distribution  which 


was  consistent  with  the  S R-X  and 


&  ZT  tests  (see  Fig.  19)  and 


at  the  same  time  fitted  the  damage  sensitivity  curve  of 


/S  6 4  >  tc  r  /  ^  c  -  Go  -  tO  J  H  s  J  [89] 

The  deflection  distribution  for  this  set  of  constants  as  well  as  other 
Sets  are  given  in  Figures  28,  29.  The  theoretical  curve  using  the  above 
mention^  constants  is  compared  with  the  experimental  results  of  Baker, 
et.  al.  in  Fig.  30.  Using  _the  same  set  of  constants  curves  of  the  non- 

dimensional  impulse  function  r/t\fjd}  as  a  function  of  V/c>  and  Pft  for 

all  values  of  were  calculated  and  plotted  in  Fig.  31,  32. 

The  effect  of  especially  for  short  cylinders  in  which  u/p<i  is  very 

large.  For  longer  cylinders  with  >  2-  the  effect  of  &/t  is  very 

small.  As  the  decreases  the  effect  of  P't  increases.  Note  that  the 
values  of  =•  tco^  /o  used  in  Fig.  32  represent  a  vast  range.  For 
P/t  >  too  there  is  no  effect  of  0/t  whatsoever. 


K  -  /  J  C  -  I ;  Co  “  /O, 


H  -  3 
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Fig. 30  Comparison  of  Theoretical  and  Experimental 
Results  for  Damage  Sensitivity  of  a 
Cylindrical  Shell  -  L/D=2,D/t=85 
x  Experimental  Points 
(Shell  is  clamped,  Rigid-Perfectly  Plastic 
Theory  is  used) 

180°  cosine  loading 
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